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Abstract 

The emergence of collective motion, also known as flocking or swarming, in groups of 
moving individuals who orient themselves using only information from their neighbors 
is a very general phenomenon that is manifested at multiple spatial and temporal scales. 
Swarms that occur in natural environments typically have to contend with spatial disorder 
such as obstacles that hinder an individual's motion or communication with neighbors. We 
study swarming particles, with both aligning and repulsive interactions, on percolated net- 
works where topological disorder is modeled by the random removal of lattice bonds. We 
find that an infinitesimal amount of disorder can completely suppress swarming for parti- 
cles that utilize only alignment interactions suggesting that alignment alone is insufficient. 
The addition of repulsive forces between particles produces a critical phase transition from 
a collectively moving swarm to a disordered gas-like state. This novel phase transition is 
entirely driven by the amount of topological disorder in the particles environment and dis- 
plays critical features that are similar to those of 2D percolation, while occurring at a value 
of disorder that is far from the percolation critical point. 

Introduction 

Collective motion of self propelled individuals is a well studied emergent phenomenon yj that 
spans many different length and time scales from biopolymers on a bed of molecular motors 0, 
swimming bacteria birds and fish [|5l[6l to people "moshing" at heavy metal concerts [|7]|. 
Within the literature that is aimed at studying collective motion in systems of self-propelled par- 
ticles [[8l[9l[lOl[IIl[l2l[l3l, a main underlying assumption has been that the environment, where 
the particles exists, is continuous, isotopic and ordered. In the natural world there are many 
examples of disordered environments where collective motion can exist. Examples include bats 
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which navigate natural caverns via echolocation, schools of fish that maneuver through dense 
kelp forests, microbial colonies that move about in heterogeneous soil lfT4l . crowds of peo- 
ple that are evacuating a building ifTSll and traffic flow in major cities lfT6l . Given that natural 
environments can be intrinsically disordered, it is interesting to consider how self propelled in- 
dividuals maintain an organized state of collective motion without knowledge of a global "road 
map". What are the necessary physical mechanisms that facilitate the collective movement 
of particles across a intrinsically disordered network? What is the role of "thermal" noise in a 
swarm that moves through a disordered environment? In this manuscript we provide insight into 
these questions for the first time, by studying a two dimensional system in which we represent 
the disordered environment as a percolated lattice that our self-propelled particles inhabit. 

Model & Simulation 

To study swarming behavior in the presence of topological disorder we implement a Monte 
Carlo lattice gas model ^ \IT^ that consists of A^^ interacting mobile particles that occupy a 
2d periodic triangular lattice with L'^ = N lattice sites. The particle density is defined as 
p = Np/N and in this model it can, in general, be larger than the density of lattice sites. The 
dynamics of our model [see supplementary information for more details] are such that, at any 
given time step, particle k moves along any of the six lattice bonds (with direction unit vectors 
u^) with a velocity v^^^ whose magnitude is a constant defined to be unity [bond length/time 
step] (Fig.[T]\). Particle moves are controlled via a standard Monte Carlo procedure with Boltz- 
mann weights determined by two interaction energies, between nearest neighbor (n.n) particles, 
in our model. The first is an alignment interaction energy, (Eqn. [T]), that makes it favor- 
able for particles to orient their velocity vectors along the direction of the average velocity of 
their n.n. The magnitude of alignment is controlled by the parameter a, which is set to unity 
and defines the relevant energy scale. The second is a mutual repulsive interaction energy, 
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(Eqn.[2]), that is proportional to the difference in the local density n{i) — n(0) seen by a particle, 
where n{i) is the number of particles at a n.n. site along a lattice direction and n(0) is the 
number of particles at the particle's current site. The magnitude of the repulsive interaction is 
controlled by the parameter e and the local repulsive energy enhances the probability of moving 
in directions with the greatest drop in density. We introduce topological disorder in our lattice 
system by varying p, the probability that a bond exists between two lattice sites as in usual bond 
percolation theory ifTSll . The quantity, 1 — therefore represents the environmental disorder 
fraction. Finally, the parameter T, which enters via the Boltzmann weighting for the Monte 
Carlo procedure effectively controls the magnitude of thermal noise in our system. We quench 
the disorder in the lattice during each realization of our simulation, allowing us to study the 
effect of topological disorder and thermal noise independently and their effect on the fidelity of 
local information that particles use to navigate the disordered environment. 

Et = -au,-J2^k (1) 

k 

El = e{n{i) - n{0)) (2) 

Results 

In general, the existence of intrinsic noise in a system of mobile interacting particles can drive 
the departure from a collectively moving ordered swarm state to a completely disordered gas- 
like state state lfT9ll . In our system, we consider a different source of disorder, environmental 
topological disorder, which is controlled by the parameter p, in addition to thermal noise. Fo- 
cusing first on the thermal noise dependence of our model without the introduction of any envi- 
ronmental disorder (p = 1) (Fig.[T^), we find that there exists a order-disorder phase transition 
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near a critical magnitude of thermal noise Fig. [T]C. This transition is reflected in the sharp 
drop of the global alignment order parameter {v) (Eq.|3]) at T^. Here {v) is the mean velocity of 



the particles in our system and can range from unity for a perfectly ordered swarm moving 
collectively along a single direction to zero when particle motions are completely uncorrelated. 



The phase transition associated with this order parameter is in good agreement with previous 
studies of the classic Vicsek model (Hllll. To understand the effects of topological disorder on 
the formation of a collective swarm we fix the magnitude of thermal noise below the critical 
value, T < Tc. In doing so we hold the system in a regime where swarming would naturally 
occur without environmental disorder. First we consider the case where particles only interact 
via a local alignment interaction (Eqn.[T]). In Fig.[2}\ we see the effect of environmental disorder 
on the formation of a collective swarm as characterized by the order parameter Eqn.[3} averaged 
over multiple disordered lattice realizations, for a system without repulsive interactions. At low 
densities of particles (p < 0.2), consistent with earlier studies [O, we find that there exists 
no collectively moving swarm ({v) < 1) that involves all particles in our system even for a 
perfectly ordered lattice. At higher densities, for low values of the disorder fraction, we find 
a completely ordered swarming state that exhibits a phase transition beyond some fraction of 
missing bonds (1 — Pq) in our system. This order-disorder transition is completely determined 
by the amount of topological disorder, which is unrelated to the classic Vicsek model transition 
that is induced by thermal noised HI [191 120)|. To understand if this transition uniquely defines 
a critical disorder fraction in the thermodynamic sense, we characterized the location of the 
critical disorder fraction (1—Pq{L)) for finite systems as a function of the system size. Finite size 
scaling analysis revealed (Fig. |2^) that, in the large system size limit (L oc), the existence 
of a critical phase transition that is governed by environmental disorder is fully suppressed, (1 — 
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where |v^| = 1 



(3) 



Pq{oo) = 0). This result suggests that, in the absence of repulsive forces between particles, any 
amount of environmental disorder destroys the ability for a system of particles to collectively 
swarm. In contrast to this behavior, we now consider the effect of adding a mutually repulsive 
interaction between neighboring particles (Eqn.[2]). In Fig.[2]C, we see that the ability of particles 
to form a collectively moving swarm is significantly enhanced compared to Fig. [2}\ even for 
moderate values of the disorder fraction. At first glance it seems that the effect of adding 
a repulsive potential between particles has changed the location of where the order-disorder 
transition (1 — p^{L)) takes place along the disorder fraction axis for finite systems. Again, 
using finite size scaling, we found that there exists a true thermodynamic phase transition from 
an ordered swarm to a disordered one (Fig. [2^) that occurs at a finite value of the topological 
disorder fraction. For a specific value of the repulsive interaction, e = 10~^, and thermal noise 
T = 10"^, this turns out to be 1 — p*(oc) = 0.13(4). Furthermore, the location of the critical 
disorder fraction moves as the repulsive interaction is changed as shown in Fig.[2p). At a fixed 
value of thermal noise, the location of the critical phase transition is governed by ratio of the 
two interaction energies, e/a {a = 1), but the critical scaling behavior near the critical point 
may in fact be universal. The universality class is usually determined by the values of the critical 
exponents that govern the scaling of various physical quantities near the critical point. We first 
extended our finite size scaling analysis to compute the critical exponents associated with this 
transition for a fixed value of e = 10"^. The correlation length which measures the length 
scale over which particle motion is correlated across the disordered lattice, diverges becoming 
comparable to the system size ((^ ^ L ^ {p^{L) — p) ^ L"^/"^) as we approach the critical 
disorder fraction transition. This allows us to compute the associated exponent u. From the fit 
to the scaling law in Fig. |2^, we find that u = 1.25(6). We also measured the susceptibility 
of the order parameter to the lattice disorder fraction by measuring the peak value of the order 
parameter fluctuations as a function of system size ^,Xv = Nal{p*) ~ Lt/^ This provided 
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an estimate for the susceptibility critical exponent, 7 = 2.70(5) (Fig. SI A). 
We also performed an independent measurement of the susceptibility exponent by examining 
how the fluctuations in the order parameter, x = L^^l{p) ^ ~ P)~^ ^ scale with disor- 

der fraction 1 — p near the critical point for different values of the repulsion magnitude. We 
found that for over three decades of the repulsive interaction magnitude that the estimate for the 



critical exponent given by Fig. [ST| \ agreed very well with estimates given by the scaling near 
p*(L) (Fig. [sT^), suggesting universal behavior with respect to repulsive strength. Scaling of 
the the order parameter, {v) ^ {p^{L) — pY with respect to p, revealed the scaling exponent (5 



(Fig. [Sip), which was also insensitive to the value of the repulsive interaction. Using the hyper 
scaling relation 2/3 + 7 = dv, we were able to infer the scaling exponent u (Fig.[STp) over the 
same repulsive interaction range and found a reasonable comparison with our estimates from fi- 
nite system size scaling (Fig.[2]C). To gain perspective on which universality class best describes 
swarming in disorder, we compared our exponents to those of two dimensional percolation [21 J 
and the standard Vicsek model [l8l[20l. We found that our system follows a universality class 
that is closer to that of percolation than the Vicsek type as shown in Fig. [ST^-D). It is inter- 
esting to note that this system is more akin to the percolation type universality class given that 
the critical disorder fraction (1 — p*(oc) = 0.13(4)) for a repulsive strength of e = 10"^ is 
vastly different than that of ordinary connectivity percolation {1 — p^ ^ 0.66)) for a triangular 
lattice [22J. This comparison suggests that swarming phenomena in such systems are extremely 
sensitive to the effects of ordinary percolation far from the percolation critical point while re- 
taining the critical behavior of a percolation type system. It is also relevant to note that this 
disorder induced phase transition occurs in a self-propelled collective swarm rather than being 
a field driven transition as is common in condensed matter systems. 

Our results show that a finite amount of repulsion enhances swarming, suggesting that there 
might be an optimal degree of repulsion depending on the system parameters. In Fig. |3}\ we 
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examine the effect that repulsion has on the order parameter in the presence of different amounts 
of disorder. In general, for finite disorder, we found that repulsion has a non-monotonic effect 
on the ability for particles to form collectively moving swarms. When the magnitude of the 
repulsion is small, we find that there is no significant enhancement of swarming ability above 
what is seen for systems without repulsion. As the repulsion is increased, we find that there is 
a maximal enhancement of the order parameter, represented by the peaks in Fig. [3}\. Beyond 
this maximum, we find that, for strong values of repulsion, there is a suppression of collec- 
tive behavior, indicating the existence of an optimal repulsive interaction for a given disorder 
fraction. 

To gain insight into this non-monotonic behavior, we look at the instantaneous disorder averaged 
mobility fraction, which measures the fraction of particles that are not temporarily stuck at a 
defect, for various values of 1 — p (Fig[3p). Increasing the repulsion magnitude, we find that 
the mobility fraction also increases for a fixed disorder fraction. Visually we can confirm that 
particle motion is becoming less hindered by the existence of lattice defects as shown in the 



simulation snapshots Fig. [S3p and D when compared to Fig. [S3) \ and C. This result is also 
consistent with the increase in the critical disorder fraction as e is increased at fixed thermal 
noise (Fig. [2p). As we approach the maximal repulsive magnitude, we find that the mobility 
fraction saturates at fixed disorder fraction (Fig.|3p) and the critical disorder fraction saturates 
as well (Fig. [2^)). Increasing the repulsive interaction past the optimum value (e ^ 10°), we 
find that the mobility fraction remains maximal for 1 — p < 0.6 (Fig. [3^), while the order 
parameter in this regime is greatly reduced (Fig. [3}\) and there exists no collective state even 
for a disorder fraction equal to zero. It is interesting to point out for a large disorder fraction 
(1 — p > 0.66), there is a slight dip in the mobility fraction for repulsion near e > 10"^. This is 
due to the geometric percolation transition for a triangular lattice, which has a critical disorder 
fraction ofl—pc = 1 — 2 sin(7r/18) 0.652 [22J. Near this transition, we find that the particles 
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become trapped in local clusters in the lattice and have a slightly higher mobility fraction for 
1 — p > I — Pc At large values of repulsive strength (e ^ 10^) the mobile fraction is greatly 
reduced at even small values of the disorder fraction, because particles are forming a gas-like 
state that interacts with defects more frequently such that the mobility scales as the number of 
missing bonds in the lattice (/x) ^ 1 — p for p < 1 as opposed to (/x) ^ 1 for 10~^ < e < 10^. 
To summarize our findings, we present the phase diagram for fixed thermal noise (T < T^) 
defined by the magnitude of the order parameter (Eqn. |3]) in Fig. [3]C. In general, we find that 
for both large values of the disorder fraction and the repulsion magnitude that swarming is 
suppressed. Moreover, for a fixed value of the disorder fraction we see that scanning through the 
repulsion magnitude (Fig. |3]C vertical axis) from low to high takes the order parameter through 
a maximum, as in Fig. |3p. At larger values of thermal noise, we find that the broad enhanced 



region is washed away (Fig. [S2p). It is interesting to note that there is a weak enhancement 
effect from thermal noise, which can be seen in Fig.[S2|\-C. We find that thermal noise (T < Tc) 
allows the system to anneal in the presence of topological disorder, which pushes out the order- 
disorder transition (even in the absence of repulsive interactions(Fig. |S3^) to slightly larger 
values of the disorder fraction (Figj2p)). Here, due to thermal noise, particles are susceptible 
to random fluctuations in their ability to align with their neighbors thus allowing the system to 
behave more ergodically. 

Conclusion 

Navigating an intrinsically disordered environment without a priori knowledge of the topology 
of the environment possess a significant challenge for individual organisms. The ability to form 
a collective intelligence in the form of a swarm can greatly reduce the navigational complexity 
in a disordered environment. A swarm utilizes physical interactions between individuals as a 
means of a primitive collective perception in order to overcome environmental obstacles. The 
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Plasmodium of Phyarum polycephalum when placed in a two dimensional maze uses a col- 
lective cellular oscillations to find the shortest path between food sources [|23l|24]l. Microbial 
colonies navigate tiny interconnected channels formed by roots of plants [145. How do these 
real systems deal with the disorder? How do they communicate local information about their 
environment to their neighbors? Topological disorder is also intrinsic in social networks, the 
internet and scientific citation networks ||25]| . One may ask the question what is the critical local 
degree of connectivity that will allow for the entire population to act holistically [|26l[27[|25ll ? 
In robotics it is becoming more common to utilize robotic drones to explore dangerous envi- 
ronments in place of humans. What are the necessary interactions that groups of robotic drones 
must possess in order to perform tasks such as search and rescue in an environment where the 
local topology may not be known a priori [|28ll ? Systems of self propelled particles using local 
nearest neighbor alignment individual particles can form a collectively moving swarm ISllll. 
However, in the presence of environmental topological disorder alignment alone is insufficient 
for particles to form a swarm. In order for these systems to navigate disorder they must possess 
information not only from local neighbors directions of travel but also the local surrounding en- 
vironment. Repulsive forces allow particles to communicate local topological features to their 
neighbors and restores the ability of the individuals to form a swarm that collectively navigates 
the intrinsically disorder environment. We have shown there exists a new type of dynamical 
phase transition driven by environmental disorder and that the ability for agents to collectively 
move in these disordered environments requires the presence of local repulsive forces between 
collectively moving particles. Collective motion in these disordered environments can be opti- 
mized by tuning the magnitude of the repulsive interaction for a given amount of disorder. In 
nature there may exist evolutionary pressures that select for better swarming ability within a 
group of individuals. It is interesting to speculate whether organisms that routinely deal with 
disordered environments have evolved mechanisms that effectively mediate repulsive interac- 
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tions similar to the one we have studied here. 



Supplementary Information 

Model-Monte Carlo 

In this section we detail our simulation methods. Lattice configurations were generated ran- 
domly for a given value of the disorder fraction 1 — p and system size N = L^. During any 
particular realization of the simulation, the bond topology was kept fixed. Particle updates were 
carried out at each time step by selecting a particle at random. The particle's direction of mo- 
tion along one of the six lattice directions (i = [1,6]) was updated by computing the probability 
to move Pi (Eqn.S [TJ in a given direction, using the local alignment field (Eqn. [T]) and re- 
pulsive particle-particle number density (Eqn. [2]). In general each of lattice sites has a number 
> particles which maybe be overlapped and can point in any of the six lattice directions. 
Topological disorder is coded into the update probability. Pi by the condition that when a bond 
exists between two neighboring lattice sites then r]i = 1 and conversely when a bond is missing 
r]i = 0. Trial moves were selected using a standard METROPOLIS Monte Carlo method where 
particle moves were accepted/rejected by comparing a randomly generated number < r < 1 
to a mapping of the move probabilities (Eqn.[T]) to the interval [0,1]. Lattice directions that have 
a missing bond are assigned a move probability equal to zero. If a particle encounters a broken 
bond it will wait there until the next move update. The total number of particle updates to reach 
steady state are system size dependent and vary from 10^ for A^^^ = 16^ to 10^ for A^^^ = 128^. 
All averaged quantities are disorder averaged over different realizations of the lattice network 
at fixed values of 1 — p. 



Pi = P{ui) = -^ry^J^exp - - arjjUi • 



X exp(-/3£;[), (1) 



10 



where Z = Y.n=i ^n- 
Model-Discussion 

To understanding how topological disorder affects the ability for particles to form a swarm we 
couple the disorder directly to the dynamics of the particles local information. This is achieved 
by allowing particles to only communicate with their n.n when a bond exists between them, 
determined by the conditional r]i. Treating disorder this way, missing bonds are represented 
as infinite barriers that forbid local information from influencing particle motion as well as 
inhibiting motion along that bond. When 1 — p > the available amount of local information 
is limited as there is an average coordination number (z = , z = 6 for a triangular lattice) 
available to each particle. In Fig. [T|\ we see an example of how a missing bond affects the 
next possible move that a particle can make. Consider the center particle in Fig. [T|\ and for 
the moment we will set e to zero. If all lattice bonds were present (i.e all bonds are black) the 
particle would determine that the most probable direction would be to continue on along its 
current direction given by the grey horizontal arrow (along the U4 direction). In the case when 
the one bond is missing (grey) the particle has limited information about the local alignment 
field and will most likely choose to move along the white arrow, which is along the U3 direction. 
Furthermore, particles can passively run into a dead end where a bond is missing from the 
lattice. Once this happens motion ceases for that particular particle until a new direction is 
computed using Eqn. [T] that takes the particle away from the broken bonds. For this reason 
broken bond directions are not allowed to be chosen by the Monte Carlo to avoid particles 
becoming permanently stuck at a broken bond. 
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Finite Size Scaling 

To provide estimates for critical exponents we measured the system order parameter, (v) and 
the fluctuations of the order parameter, a^, for various system sizes at a fixed density, p = 1. 
Using the finite size scahng functions below we extracted both the critical disorder fractions 
threshold and critical exponents (Fig. [2^). To locate the critical disorder fraction in Fig. [2^ 
we use the fact that the correlation length will scale as the system size near the critical disorder 
fraction, ^ L, 

p*(L)-p*(oc)^L-^ (2) 

where X = 

The susceptibility is related to the fluctuations of the order parameter near the critical disorder 
fraction, 

X.= j^[{v')-{vr]=L'al (3) 
The peak of the susceptibility also scales as a function of system size (Fig. [ST|\) , 

X„(p*(L),L)^L^/^ (4) 

To check the universaUty of these exponents over the range of repulsion magnitudes we fit both 
the order parameter and the susceptibility near the critical disorder fraction, p*(L)(Fig. [sT^ 
and C). 

{v)^{p-p*{L)f (5) 

and 

Xv-ip-p*{L))-''^L'al (6) 
Using these estimates we used the hyper scaling relation, 

2/3 + 7 = o?i/ (7) 

to compute u when varying the repulsive interaction(Fig. [STp). 
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Figure 1: (A) Diagram of a lattice unit cell with particles (red) occupying lattice sites. Solid 
black bars are occupied lattice bonds and grey bars are deleted bonds. The blue triangles rep- 
resent current particle velocity directions. White/Grey arrows indicate some potential lattice 
directions along which the center particle can move at the next time step. In the absence of 
any repulsion, if the deleted bond (grey bond) were still present, the updated direction of the 
center particle would most probably be along the grey arrow, since that is the local average 
velocity direction. However, since the bond is not present, the most probable direction is along 
the white arrow. Note that the absence of the bond not only disallows motion in that direction 
but also inhibits transfer of information, (see supplementary information)(B) A snapshot of a 
typical simulation showing finite sized groups of particles which are collectively moving in a 
disorder free (1 — p = 0) lattice. (C) Order-disorder transition driven by the thermal noise in 
a system of particles (p = 1.0 and 1 — p = 0) without repulsive interactions for a system size 
of = 1024. The blue dash line indicates the location of the critical transition temperature, 
Tc 0^ 0.25, which was determined numerically. (D) A snapshot of a simulation of a system with 
disorder (1 — p = 0.05), where missing bonds are not drawn. Black filled circles are particles 
that are temporarily stuck at a lattice bond defect. 
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Figure 2: (A) Global alignment, in the absence of repulsion, as measured by Eqn.|3]as a func- 
tion of the disorder fraction in the lattice for different densities [p = 0.2 (black), p = 0.8 (red), 
p = 1.0 (blue), p = 2.0 (green) ]. (B) Finite system size critical disorder fraction (1 —p^{L)) for 
systems without repulsion [6 = (green filled circles)] and with repulsive interactions [e > 
(blue filled circles)]. Fitting to finite size scaling, — p*(oc)) ^ L"^/"^, predicted a critical 

disorder fraction for systems without repulsive interactions (black dashed line) that is zero and 
for systems with repulsion (red dashed line) is found to be 1 — p*(oc) = 0.13(4)(black solid 
line). The critical scaling exponent u in the case for repulsive interactions is 1.25(6). (C) Global 
alignment for a system of particles that interact with both a local alignment and repulsive fields 
for p = 1.0 and e = [5.0 * 10"^ (black), 10"^ (red), 5.0 * 10"^ (blue), 10"^ (green)]. The lo- 
cation of the order-disorder transition is pushed out allowing particles to swarm in the presence 
of moderate environmental disorder. (D) The location of the critical disorder as a function of 
temperature at various values of the repulsive interaction strength (see inset). For low values 
of repulsion the critical disorder changes very little over temperature. At intermediate repul- 
sion values we see there is an enhancement effect with increasing thermal noise. At repulsive 
strengths near maximal (Fig. [Tp) we find litti^ change of the critical disorder fractions over 
entire range in T. All data shown in (A), (C) and (D) is for a system size = 1024. All data in 
(A), (B) and (C) presented here was obtained at a thermal noise of T = 10~^ 




0.1 0.2 0.3 0.4 0.5 

Disorder (1 - p) 

Figure 3: (A) Behavior of the order parameter for particular values of the disorder fraction 
[p = 1.0 (red), 0.90 (blue), 0.80 (green)] clearly shows a non-monotonic dependence on the 
magnitude of the order parameter (log horizontal axis) with a central maximum for curves 
with p < 1. (B)The ensemble averaged instantaneous mobility of particles as function of 
environmental disorder. For smaller values of the disorder fraction, the effect of increasing 
the repulsive energy significantly increases particle mobility allowing particles to swarm. As 
the disorder fraction is increased, the magnitude of the mobility generally decreases as expected. 
Interestingly, for values of the disorder fraction near the connectivity bond percolation threshold 
(1 — Pc ^ 0.66) we find a significant dip in the mobility of particles at large values of the 
repulsive energy parameter (e > 10^). This signifies that particles are becoming trapped in 
local disconnected clusters in the lattice. (C) Phase diagram defined by the magnitude (color 
bar) of the order parameter (Eqn. [s]) for the two parameters e and 1 — p dXT = 10"^. The 
large peak centered near log(e) ^ — 1 shows the effect of repulsion has on the persistence of 
swarming the presence of disorder. 

17 



A 



10 



^10' 

bJO 
O 

^ 10^ 



10" 



7 = 2.70(5) 



log(L) 




0.1 



^ ^ -r ¥ * 

■ ■ ■ ■ ■ ■ ■ ■ ■^F ■ ■■■■■■■ 



0.05 



-3 -2.5 -2 -1.5 -1 

log(e) 




-3 -2.5 . -2. .-1.5 




-3 -2.5 -2 -1.5 

log(e) 



Figure S 1 : (A) Finite size scaling fit (red dashed line) of the susceptibility (blue filled circles) 
at finite repulsion (e = 10~^) as a function of system size. The critical exponent, 7, associ- 
ated with the scaling of the susceptibility is shown. (B) Extracted critical exponents with error 
bars for 7 (blue circles and black bars) from scaling analysis near the critical disorder fraction 
transition (1 — p*(I/)) plotted against the log of the repulsion magnitude. The grey dashed line 
is the average of the these values while the red line is the predicted exponent for percolation 
and the green line is the predicted exponent for the standard Vicsek model. (C) Filled circles 
are the extracted exponents, /3 for the scaling of the order parameter near the critical disorder 
fraction [see supplementary material]. The dashed lines have the same color arrangement as the 
previous figure in the series. (D) Filled circles are the estimates for the critical exponent, z/, as- 
sociated with the correlation length, plotted over the magnitude of the repulsive interaction [see 
supplementary material]. Dashed lines indicate the averaged exponent and expected exponents 
from percolation and the standard Vicsek model. The color arrangement is the same as in the 
previous figures in this series. 
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Figure S2: (A) Increasing thermal noise, for system without repulsion, but below the critical 
value there is an annealing effect that pushes the order-disorder transition to slight higher values 
of disorder fraction. (B)The combination of both the repulsive energy and the temperatures near 
Tc (T = 10"^) both help the ability for swarming behavior near the critical bond occupation 
probability for e/a <C 10"^ as compared to Fig.|3}\. (C) A phase diagram for larger thermal 
noise magnitude (T = 10"^) near the critical value Tc (compare to Fig. [T]C). The prominent 
peak near log(e) ^ — 1 has diminished significantly, but there also is a slight thermal noise 
effect which restores swarming for weak values of the repulsion magnitude (log(e) ^ —3) and 
small values of the disorder fraction. 
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Figure S3: Snapshots of a simulations with and without repulsive energy. (A) In early simula- 
tion times without repulsion, particles scatter off locations where bonds are missing and begin 
to form compact single file lines. (B) With repulsion (e = 10~^) in early times, particles form 
extended groups that move together around defect bonds and fill in gaps caused by defects. (C) 
At late times without repulsion, particles form ^8 extended single line that finds the a defect free 
lattice direction. Black particles are stuck at a broken bond indefinitely in this case. (D) In late 
times with repulsion, we find an collectively moving ordered swarm which can avoid defects 
and remain ordered by moving around particles which are temporarily stuck at a broken bond. 



